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Abstract 

The operator product expansion (OPE) for heavy-light-quark pseudoscalar mesons (D- 
mesons and B-mesons) in medium is determined, both for a moving meson with respect 
to the surrounding medium as well as for a meson at rest. First of all, the OPE is given 
in terms of normal-ordered operators up to mass dimension 5, and the mass of the heavy- 
quark and the mass of the light-quark are kept finite. The Wilson coefficients of such 
an expansion are infrared (IR) divergent in the limit of a vanishing light-quark mass. 
A consistent separation of scales necessitates an OPE in terms of non-normal-ordered 
operators, which implies operator mixing, where the IR-divergences are absorbed into the 
operators. It is shown that the Wilson coefficients of such an expansion are IR-stable, 
and the limit of a vanishing light-quark mass is perfomed. Details of the major steps for 
the calculation of the Wilson coefficients are presented. By a comparison with previous 
results obtained by other theoretical groups we have found serious disagreements. 
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I. INTRODUCTION 



The discovery of the charm-quark in 1974 by detecting the J/\l/-particle l|, |2| has 
completed the second quark-generation and was a manifest triumph of the quark- 
model. While J/\l/ is a bound state of a charm-quark and an anti-charm quark, 
J/\I' = cc, so-called open-charmed mesons were discovered soon afterwards: = 
cd, = cu, and their corresponding anti-mesons = dc,l^ = uc. The open 
charmed mesons are much easier accessible experimentally, because the D-mesons 
are the lightest particles which contain a charm-quark and their lifetime is a few 
orders of magnitude larger than J/^>: the mass of J/\E' is 3096.9 MeV and the 
lifetime is 7.06 x 10~^^ sec, while. 



and the lifetime 1.04 x 10 sec 



'or instance, the mass of D^-meson is 1869.6 MeV 



|3|. 

Whereas charmed mesons in vacuum were studied thoroughly ever since, noth- 
ing is known so far from experimental side about the properties of these mesons 
in medium. Nowadays, experiments are initiated to study charmed mesons embed- 
ded in a hadronic medium. In fact, the upcoming accelerator FAIR (Facility for 
Antiproton and Ion Research) at GSI (Gesellschaft fiir Schwer-Ionenforschung) in 
Darmstadt /Germany ^| offers the opportunity to study charmed mesons in dense 
nuclear matter. Especially, the CBM (Compressed Baryonic Matter) collaboration 
intends to study the near-threshold production of D-mesons and J /'^ in heavy-ion 
collisions, while the PANDA (anti-Proton ANihilation at DArmstadt) collabora- 
tion will focus on charm spectroscopy as well as on charmed mesons produced by 
anti-proton annihilation in nuclei. 

This growing interest can also be motivated from in-medium modifications of K- 
mesons, that is an expected down-shift of K~ and an up-shift of with increasing 
density [5i]. There is a similarity between D-mesons and K- mesons in respect to 
their quark structure: K~ = su corresponds to , while = us corresponds 
to . Thus, we conclude the qualitative fact, that both strangeness and charm 
might be regarded as probes of the in-medium situation. Moreover, the expected 
in-medium modifications_of D-mesons might have a considerable impact on normal 
J/\I' suppression, e.g. [6], and open-charm enhancement, e.g. [jl, in heavy-ion 
collisions. 
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Despite the upcomming experiments CBM and PANDA, theoretical investiga- 
tions of charmed mesons, both in vacuum and in medium, are fairly rare. One rea- 
son is, that the well-established perturbation theory of Quantum Chromodynamics 
(QCD), the fundamental theory of strong interactions, cannot be applied, because 
ihe momentum transfer Q among bound quarks is small Q ~ Aqcd — 245 MeV 
8|. Thus, the running coupling constant of QCD as{Q) becomes large, and non- 
perturbative approaches are compelling in order to derive the physical properties of 
mesons from fundamental principles. 

Especially, the infrared regime of QCD is complicated due to the still unre- 
solved phenomenon of confinement. One aspect of thereof is the non- vanishing of 
normal-ordered products of quark and gluon field operators taken between the non- 
perturbative, so-called physical, groundstate |vac) of QCD, e.g. the chiral quark- 
and gluon-condensate: 

(vac| : g^'g'" : |vac) ^ 0, (1.1) 



(vac 



|vac)^0. (1.2) 



Here, : O : denotes normal-ordering of a composite operator (9, is the quark 
field operator of color index a = 1,2,3, and is the gluon field strengh tensor, 
where /i, = 0,1,2,3 are Lorentz-indices and A = 1,...,8 is the Gell-Mann index; 
troughout the article the Einstein convention is used. Vacuum expectation values 
like (11. II) and (II. 2p are called vacuum-condensates, a term which refers to the fact 
that the physical ground state of QCD is not "empty" but contains colorless and 
chargeless states of strongly bounded (condensed) quarks and gluons. In contrast, 
in the perturbative vacuum |0) of QCD the normal-ordered condensates vanish, e.g.: 
(0| -.q^q": |0) = and (0| : G^^G^^"" : |0) = 0. 

At finite baryonic densities the vacuum-condensates are generalized to in- 
medium-condensates, where the in-medium state is denoted by i.e. 

{n\:rq':\^), (1-3) 
{n\:G^^G^^'':\n). (1.4) 

The ground state of a hadronic medium can be approximated as a Fermi-gas of 
nucleons: \Q) = |vac) + |N(/c)). It is a formidable challenge to evaluate vacuum 



and in- medium condensates from first principles of QCD. So far, approximative 
solutions have been obtained (e.g. lattice gauge theory, Wilson loop expansion, 
instantons, etc.), but a comprehensive understanding of the confinement problem is 
far out of reach. Therefore, the question arises how fundamental QCD parameters 
may be related to the vacuum and in-medium properties of mesons. 

One of the best methods among the non-perturbative approaches are the QCD 
Sum Rules (QSR), which link hadronic observables like mass, decay- width and 
coupling-constant of the hadron under consideration to fundament al p arameters 



of QCD. The QSR have first been developed for mesons in vacuum 



9l-|ll[| . and later 

been generalized to the case of mesons in medium [I2l4l5l|. see also 16j for the case 
of finite density and temperature. The decisive advantage of this approach is, that 
QSR consider the existence of vacuum-condensates and in-medium-condensates as 
a given fact, instead of determining their numerical values by first principles of QCD. 
In this way, QSR are not concerend with the problem of confinement, but assume 
that confinement exists and they are neither concerned with the involved structure 
of the groundstate |vac) nor the in-medium state 

Thus, once a very few parameters, i.e. numerical values of condensates, are 
given, physical properties of a large variety of hadrons can be predicted. The link 
between hadronic observables and fundamental parameters of QCD is given by a 
dispersion relation of n-point correlators which are non-local products of n opera- 
tors: Oi{xi)02{x2)...On{xn)- Here, we consider so-called current- current correlators 
j2{y), which are two-point correlators. 

The analyticity of the current- current operator allows to relate the time-like re- 
gion (hadronic part) with the space-like region (QCD part) of that operator. While 
the hadronic side of the dispersion relation is parametrized by means of hadronic 



observables, the Wilson Operator Product Expansion (OPE) [171] is applied on the 



QCD side. In general, the product ji{x)j2{y) diverges in the limit x — )■ ?/, e.g. 
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43| . The Wilson OPE allows to decompose the product of non-local opera- 
tors in terms of a series of regular local operators O and so-called Wilson coefficients 
Ck which are divergent for x — )■ ?/, i.e. 

Mx)j2iy) = fyu (^) O, . (1.5) 



The point of the OPE (11. 5p is that the fields of the currents are separated into a 
hard and a soft part. The hard part is proportional to the unit operator and can 
be treated perturbatively, i.e. can be evaluated with respect to the perturbative 
vacuum |0). The spatial dependence of the soft part can be Taylor expanded and 
leads to local condensates in the medium state 

In general, the currents of the correlator are interpolating fields, which carry 
the symmetries of the specific hadron: spin, isospin, parity, charge and the valence 
quark content. More specific, in case of D-mesons or B-mesons we need to analyze 
correlators of heavy-light-quark currents. OPE's of such heavy-light-quark current- 
current correlators have been evaluated since the early days of QCD sum rules, see 



pioneering investigations 



20|-l27|. 



One of the main peculiarities within the calculation of the OPE for a pseudoscalar 
heavy-light-quark meson correlator is the absorption of infrared mass divergences 
which occur in the Wilson coefficients. In order to render the OPE finite, these 
divergences have to be absorbed into the condensates by introducing non-normal- 
ordered condensates, like 

(nrgin), (1.6) 

{Q\G^,G^^-\Q), (1.7) 

instead of normal-ordered ones. The unique mathematical scheme behind is the 
approach of operator mixing which allows for a consistent separation of scales. This 



problem has been worked out in detail in 27|, |28|, and accounts for a perturbative 



piece of the condensates, cf. |25|, l27H29| for a detailed discussion of the vacuum case. 

At non-vanishing baryonic densities the OPE differs from the vacuum case and 
additional mass singularities occur. In order to treat them consistently one has to 
find additional expressions which reproduce the vacuum limit but also render the in- 
medium OPE infrared stable. A crucial point is the mixing of different condensates 
under this procedure. The in-medium OPE for open charmed mesons has been 



evaluated in 
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32|. Unfortunately, the results and Wilson coefficients presented 
in these references significantly differ from each other. Moreover, it is not clear in 
which way the authors have dealt with the infrared mass singularities. However, the 
knowledge of the correct OPE is compelling for a reliable prediction of in-medium 



properties of D-mesons within the QSR approach. Recently, in [33|, [SJ] an OPE 
for heavy-hght currents in medium has been presented where the operator mixing 
and cancellation of infrared simgularities were correctly taken into account, but no 
further details of the involved evaluation of the OPE were given. In view of the 
progressing in-medium D-meson physics, both experimentally and theoretically, it 
is timely to present a transparent and thorough calculation of the OPE for D-mesons 
in matter and a comprehensive representation of the techniques which have to be 
applied. 

The pseudoscalar B-mesons are also heavy-light-quark systems and their quark 
structure is: = ub, = bu, B^ = db and B^ = bd. So far, the only application 
of the QSR approach has been performed in 33| in order to determine in-medium 
modifications of B-mesons. Since the formalism of OPE and QSR can easily be 
extended from D-mesons to the case of B-mesons, we will incorporate these mesons 
in our investigation. 

The paper is organized as follows: The OPE in terms of normal-ordered con- 
densates is calculated in Section [Tll Non- normal-ordered condensates and, therewith 
associated, the operator mixing and the absorption of infrared divergences is dis- 
cussed in Section UTTl The IR-limit, i.e. the limit of a vanishing light-quark mass, is 
considered in Section IIVI In Section |V] the OPE is given for the physical situation 
where the meson is at rest with respect to the surrounding medium. Furthermore, 
a Borel transformation is performed. A comparison with the results given in several 
publications is given in Section IVII The summary can be found in Section I VI II 



II. OPE IN TERMS OF NORMAL-ORDERED OPERATORS 
A. OPE with IR-divergent Wilson coefficients 

Let us introduce the current-current correlator for D-mesons, which is defined as 
Fourier transformation of the expectation value of the time-ordered products of two 
currents: 

n(g) = t J rf^a:e'^^(fi| Tj(x)jt(0) , (2.1) 
where T denotes the Wick time-ordering, = (go, q) is the four-momentum of the 
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D- meson, and is the groundstate of hadronic matter. Application of Wick's 
theorem to the correlator (12. ip naturally yields the OPE in terms of normal-ordered 
operators and can be written as 

n(g) = CoI+J2'^^ {n\ : : \n)^^\ (2.2) 

i 

where X is the unit operator, the notation means tree-level matrix elements, and 
the bar denotes Wilson coefficients which correspond to tree-level matrix elements. 



cf. the OPE of heavy-light-quark currents in vacuum [28|. From now on we drop 
the label i.e. normal-ordered condensates are always on tree- level throughout 
the work. 

As mentioned in the Introduction, the currents in fl2.2l) are constructed such that 
they contain the quantum numbers of the particle under consideration, i.e. spin, 
isospin, parity, charge and the valence quark content. Accordingly, for D- mesons 
they are given by 

jo+ix) = : id75C : , (2.3) 
jD-ix) = -Ac^^d:, (2.4) 
jz)o(x) = : iu'j^c: , (2.5) 
j_o(a;) = : ic75U : . (2.6) 

Due to = jo- and jj^o = jjf "we obtain, by means of j{x) = 

exp{iPx) j{0) exp{—iPx) where P is the momentum operator, the relations 
^D+{(l) = n£)-(— g) and Ilr)o{q) = n-^o(— g), respectively. Furthermore, in isospin 
symmetric nuclear matter a replacement u ^ d does not change these four corre- 
lators. Thus, it is sufficient to consider the current-current correlator (12. ip with 
the current operator of the meson given by Eq. (12. 3p . since the other three 
correlators Hd-, and H-^o, can easily be deduced from IId+- 
Similarly, for the B-mesons the currents read: 

jB+ix) = -.ib-f^u:, (2.7) 

jB-ix) = : iu'j5b : , (2.8) 

j^o(x) = ■Ab'j^d:, (2.9) 

j_o(x) = -.id^^b: . (2.10) 



Obviously, we have Jb- = and j-^o = j^,, and IlB+{q) = IIb-{—(i) and IlBo{q) = 
n-go(— g). Therefore, the OPE for the B-mesons can be deduced from the OPE of 
B^. Furthermore, since the OPE of can be obtained from the OPE of by the 
replacements u ^ d and 6 — )■ c, we conclude that the OPE of all heavy-light-quark 
pseudoscalar currents can be deduced from the OPE of the meson; therefore, 
from now on we will drop the explicit notation . 

In this Section, we will determine the Wilson coefficients of (12. 2p . The OPE 
contains all operators up to mass dimension 5: 

mass dimension , (2-11) 

mass dimension 3 , (2.12) 

mass dimension 4 , (2-13) 

mass dimension 4 , (2-14) 

mass dimension 5 , (2.15) 

mass dimension 5 , (2-16) 

where the flavor of the quark-flelds are either charm or down quarks, the Dirac 
indices are denoted by i,j = 1,2,3,4, and the operators have to carry a colorless 
structure; for notation see Appendix |Al In what follows, we will calculate all Wilson 
coefficients at first non-trivial order in powers of the QCD coupling constant, except 
for the unity operator where we include corrections. 

As we will see in this Section, the Wilson coefficients of (12.21) are divergent in the 
limit of a vanishing light-quark mass — 0, which is a so-called IR-divergency. The 
reason for this divergence is that a proper factorization of short and long distance 
contributions in the OPE requires the calculation of matrix elements at the same 
order as the Wilson coefficients. This special issue will be the topic of the next 
Section. In this Section we will consider the OPE (12. 2p . i.e. in terms of normal- 
ordered operators. 

By inserting the current (12.31) into (12.11) and applying Wick's theorem to the 



X 



cd ~r)da b 



ttC^mi' r'ca b 
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current-current correlation function n(g) we obtain: 

n(g) = nW(g) + nf (g) + n(2)(g) + n(^)(g), 



(2.17) 



n(°)(g) = 


-I j rf^xe*''^(fi| :Trc,D (75 ^d(0, x) 75 ^e(^ 


,o)):|n), 


(2.18) 


n?^(g) = 


y rf^xe^''"(fi| :5(x)75^c(x,0)75C^(0) : \n) 


; 


(2.19) 


= 


j rf^xe^''^(^]| : c(0)75 5d(0,a;)75c(x) : 




(2.20) 


n(^)(g) = 






(2.21) 



where the notation Trc,D nieans trace over Dirac- and color-indices. 

The point of the whole OPE is that the quark and gluon fields are separated into 
a hard and a soft part. The hard part can be treated perturbatively. The soft part 
cannot be calculated in this way, but being soft one can Taylor expand the spatial 
dependence of the fields and relate it to local condensates. 

Accordingly, the term n^°) is decomposed into a perturbative part, which is pro- 
portional to unit matrix and is treated by standard perturbation theory of QCD, 
and a gluonic part where the gluon fields are soft: 

nW(g) = nP^'^*(g) + ng(g) . (2.22) 

The terms n[^], H^^'*, Hc^^ and n^^-* describe the non-perturbative part of the correla- 
tor. The labels and denote the number of non-contracted quark fields, i.e. 
the number of quarks which participate in the formation of a condensate. The term 

n(4) 

concerns only the soft part of all operators. Thus there is no flow of hard mo- 
mentum. The result of the integral in ( 12.2ip is proportional to Dirac's delta-function 
5{q^) and derivatives thereof, i.e. n^^-* vanishes except for = 0. Therefore at large 
|g^| (OPE) there is no contribution from this term. 

The quark propagator in a weak gluonic background field in coordinate space 
reads (n is either c or d) 

iSn{x,y) = iS^^\x-y) 
00 „ 

+ d^^i . . . d^Zk tS^^\x - z^)gsA{zi)tS^^\zi -Z2)... gsA{z^)iS^^\zk - y) , 

k=i 

(2.23) 
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where S^'^\x — y) denotes the free propagator and 'hat' a contraction with Dirac 
matrices A = 'J^A'^. For our investigation we need the quark propagator up to order 
k = 2, given exphcitly by Eqs. (IBT3|) - (iBTSj) . 

In what follows we will determine the Wilson coefficients of f l2.18p - fl2.20p . while 
(I2.2ip does not contribute at large 



B. The perturbative part W 



ert 



The perturbative part is the Wilson coefficient of the unit operator (12. lip , and 
is given by 



(2 



Trc,D (75 SjP\p) 75 Sf\p -q))+0 (as) . (2.24) 



Note, there is neither a normal-ordering nor an expectation value, since in perturba- 
tive QCD the propagator is a usual c-number. The first term in (I2.24p corresponds 
to a one-loop Feynman diagram; see the left diagram in Fig. [H The terms of order 
O (as) are not given explicitly and correspond to two-loop Feynman diagrams, see 
the middle and the right diagram in Fig. [1] Their mathematical structure is almost 
identical to (12.280 - (12.300 . but the two gluon fields are not soft anymore but con- 



tracted to a free gluon-field propagator 



found in standard text books, e.g. 18|, |43| 



'or equal quark mass these terms can be 



The explicit solution of (I2.24p is rather involved and is usually determined by means 
of the dispersion relation 



Co = np--*(g) 



1 



ds 



ImnP^''*(s) 
s — 



(2.25) 



(mc+mj)2 

The imaginary part of the perturbative term can be determined by means of the 



Cutkosky cutting rules [45[ (for a didactical representation of the cutting rules in- 



cluding some examples see ^|), and has been calculated at the 



very first time in 



231], and later on several times by different authors, especially 20|, |26|; note that an 

23] has been corrected by [26]. In the MS-scheme 4^], 



error in the calculation of 



the perturbative part is IR-convergent [48|]. In accordance with the whole approach, 
we present our result for = in Feynman gauge, which agrees with the findings 
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FIG. 1: Feynman diagrams for the perturbative part IP'^^^{q) in Eq. (I2.24P : the fourth 
diagram is not shown but can be obtained from the right one by the replacement d -f-T- c. 
Solid lines depict the free quark propagator, curly lines stand for contracted gluon fields, 
i.e. free gluon propagator, and dashed lines denote the D-meson. The left diagram is 
for the zeroth order 0{a^), while the other two diagrams are the first order 0{as) in 
perturbation theory of QCD. 



of |20|,l26j: 



ImnP"''*(s) 



m 



2n2 



TT 
2n2 



2 7r2 



9 ■ f^c \ 



In 



3 , 
+ - In 
2 \ s 



+ In 



mt 



+ 



m; , IS 
— In 



In 



mt 



mt 



(2.26) 



Here, Li2(x) = — J dt ln(l — t) is the Spence function. We note, that the 



perturbative part must be symmetric in exchanging c-quark and d-quark. However, 
since we have presented the perturbative part in the IR-hmit m^ = 0, the given 
solution (12.261) is not symmetric anymore. 

We will not finish this paragraph without a final remark about the dispersion rela- 
tion (I2.25p . According to (I2.26p . the dispersion relation (I2.25P gives an infinite result 
in ultraviolet regions of integration, i.e. one might prefer a twice-subtracted disper- 
sion relation instead of (12.250 . Such a subtraction scheme yields polynomials in the 
external momentum q. However, after a Borel transformation has been performed, 
all polynomials disappear, that means a Borel transformation of a twice-subtracted 
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9 rud g '^'^ g 

"C^y~ "{^^^" --(^^^- 

g ^ rric ^ rric 



FIG. 2: Feynman diagrams corresponding to II'^ ' , H'^ ' and H*^ ' in Eqs. ()2.28p . ()2.29p 
and (|2.3U|) from left to right. Solid lines depict the free quark propagator, dashed lines 
denote the D-meson, and curly lines stand for the soft gluon fields; the crosses symbolize 
the creation or annihilation of a gluon by virtual particles. 

dispersion relation is identical to a Borel transformation of fl2.25p . 



C. The gluonic part II^j 

Let us turn to the gluonic part, Eqs. f l2.18p and f l2.22p . which is the result of in- 
serting the next-to-leading order propagator with the lowest order term of Eq. ( 1B12I) , 
and accounts for pure gluon condensates f l2.14p . The expressions up to order O (as) 
are given by: 

=n^^'^(g) + n^^'^(g) + n^^'^(g), (2.27) 



n^''"(g) = -'j ^M-- Trc,D (75 S^^\p) 75 S'^\p - g) ) : , (2.28) 

n^'"(?) = -^j ^(^1 :Trc,D(75^f (P)75^f (p-g)) : 1^^) , (2.29) 

n^""(?) = / (^h Trc,D (75 S^\p + q) 75 Sf{p) (2.30) 

where the quark propagators are given by Eqs. ( 1B13[) - ( IBlSp . These expressions 
correspond to the Feynman diagrams shown in Fig. [21 

First we note, that by inserting the expressions (IBlSp - (IBISP into fl2.28p - ( I2.30p 
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we recognize that these three integrals are ultraviolet finite. We also note, that 

(3) 

inserting higher orders of the quark propagator, e.g. S*^ , or the gluon fields would 
lead to higher dimensional gluon condensates, such as (: :), or, by usage of the 
equations of motion, to quark and mixed quark-gluon condensates. These terms are 
either of higher mass dimension or of higher orders in ag, so we omit such terms. 

At this stage of our investigation we would like to mention that the calculation 
of the gluonic part 11^] is performed at the one-loop level and leads to cumbersome 
expressions. Some details of the evaluation are therefore shifted to Appendix [Cl 
where we have described in detail the techniques by means of which we obtain the 
following final result for the terms fl2.28p - f l2.30p : 

1 1 



TT 




6 (g^ — m^)^ 9 (g^ — m^)^ 



(2.31) 



Here, = G^^G^^'^, and = (l,v) is the four-velocity of the surrounding 
medium. The first line in f l2.3ip is the scalar contribution which does not van- 
ish in vacuum, while the other term is a medium-specific condensate and vanishes 
at zero density. One immediately observes an IR-divergergent term oc m^^ known 
from the vacuum OPE of D-mesons, while in medium there is an additional loga- 
rithmic IR- singularity oc Inm^. These both IR singularities appear because in the 
corresponding diagram (right diagram in Fig. [2]) there are three light-quark propa- 
gators with the very same momentum, while in the other diagrams there are only 
one or two. 

We note the symmetry of f l2.28p - fl2.30p in exchanging the charm and down quark. 
However, since we have to perform the IR-limit — )■ after operator mixing, we 
have taken into account the leading IR-divergent terms only, that means m^^ and 
oc Inm^. Hence, the solution fl2.3ip is no longer symmetric in exchanging down 
and charm quark. In Appendix [C] further details are given about how the needed 
expansions in terms of a small d-quark mass destroy this kind of symmetry. 
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Furthermore, terms of the form In and In occur simultaneously and cannot 
be made small at the same time for —q^ ^ m^. However, as has been found in 48| . 
they are remnants of the large distance behavior, i.e. they originate from the small 
momentum contribution to the loop integrals. Their occurrence breaks the neat 
separation of scales, which is a necessary feature of every OPE, and must therefore 
be absorbed into the condensates. In 48|] the author argues that these logarithms do 



not occur when the Wilson coefficients are calculated within a minimal subtraction 
scheme. Moreover, they can be absorbed into the condensates if one reexpresses 
normal-ordered condensates, which naturally emerge if one applies Wick's theorem 
to fl2.ip . by so-called non- normal-ordered ones. This procedure has been known for 



a long time in vacuum, cf. 



21 



28 



29 



49| and references therein, although an 



explicit formula could not be found by us. In Section IIIII we introduce a formula 
which relates normal-ordered and non-normal-ordered condensates in matter. 



D. The term H?) 

Let us now consider the term fl2.19p . which leads to condensates which contain 
down-quarks. After expanding the light-quark fields and performing the Fourier 
transformation we obtain the expression 

fc=0 

(2.32) 

The quark fields and their covariant derivatives have to be calculated at the origin, 
i.e. X = 0. From here we can go to higher quark field derivatives or to higher orders 
in the perturbative propagator or to higher orders in the gluon field, which enters 
through the perturbative quark propagator. The quark fields are of mass dimension 
3/2. Each covariant derivative and the gluon fields enlarge the mass dimension 
by one unit. Thus working in lowest order of the gluon fields the following terms 
have to be considered up to mass dimension 5: 

n?^(g) = n?)^(g) + ni%(g) + uf^M + ni%(g) , (2.33) 
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FIG. 3: The first three Feynman diagrams corresponding to 11^ in Eq. (j2.32p . Solid lines 
depict the free quark propagator, dashed lines denote the D-meson, and curly lines stand 
for soft gluon fields; the crosses symbolize the creation or annihilation of either a gluon or a 
quark by virtual particles. The left diagram corresponds to the quark condensate (: dd :), 
while the middle and right diagrams correspond to the mixed quark-gluon condensate 
(: dgaQd :). 

where the individual contributions are 

niS(^) = (^1 ■■ < (75^f (9)75)r, d] : \n) , (2.34) 



nl%(g) = -^(^^1 : (rf^K)" (759^^f (g)75)>; : \n) , (2.35) 
= -^(^1 ■ (^^K^.)" {l^d^d'^Sj:'\qh,)ld'; : , (2.36) 

ni%(g) = (^^1 : < (755f )(g)75)!i4 : 1^^) , (2.37) 

where we note that the term fl2.35p is not present in vacuum. In Fig. [3] three 
corresponding Feynman diagrams are shown, but there is actually no one-to-one 
correspondence for the expressions fl2.34p - (12.371) and Feynman diagrams. It be- 
comes obvious, that the Wilson coefficients can be obtained on tree-level, in contrast 
to Wilson coefficients for gluon condensate which are obtained on one-loop level. 
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The calculations on tree-level are straightforward, so we just present the result: 

(g) = {n\ : dd : - : dt^.d : ll^)^^^^^ 

-{^\ : h^d : m^^^ + : ^Ta ^ : |^^) fr^^^ " 

— mi \ [q^ — m^)^ q^ — n. 

+ : d^,l5,l5.d : ^(^^^^ j^^^, ) 

-ml - 1 e"'"'^'^ g 

-.dgsagd: (^2_^2)2 + (^1 • ddsl^laQ^^.d : |^ . 

(2.38) 

At this point of the evaluation we have to keep terms oc (such terms appear when 
applying the equation of motion in some coondensates, e.g. in fourth condensate the 
term proportional to metric tensor g^^) because they will be necessary to cancel the 
infrared divergences or give finite contributions to the correlator. The limit m^ — > 
will be taken in the overnext Section after absorption of these infrared divergences. 



E. The term Ilf\q) 

Now we consider the term (I2.20p which leads to charmed condensates. In order to 
show the similarity with the expression 11^ (g) obtained in the previous Subsection, 
we will keep the light-quark mass finite ^ 0. The calculation of lie (g) implies 
a tree-level evaluation, and we just present the result: 

n(^)(g) = m■.cc■^n)^^^-{n\■.c^^,c■.\n)^^^^ 

+ m : C7.C : 1^)^^ " (^1 : cl^^^ .c : |^]) {j^^, ^ 



- m^)2 g2 - ml 



/oi -l^ l^ \n\ f ^^9'"" Amdq^'q''\ 

md 1 _ 1 e°'^'^'^ q 

-{n\ ■.cgscjgc: _ - {n\ ■. cgsi^iaG,uC ■. _ 



d) 

(2.39) 
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Obviously, there are no IR-divergent terms due to charmed condensates. In 
our study we will neglect all charmed condensates f l2.39p . This assumption can 
be justified as follows. According to the approximation of the in-medium state 
\Q) = |vac) + Ylk '^^'^ split term fl2.39p into a vacuum part and a 

nucleon part: ^ = ni'^™^ + H^^^. 

At first, let us consider the vacuum part: U^c^^^^. In vacuum, the only non- 
vanishing term in the limit m^; — )■ reads: 

n(2)vac^^) = -(vac| : c^x^^,c : |vac) (^^^ - . (2.40) 

We insert the operator mixing relation (13. 7p for the vacuum situation, and obtain 
the expression IIc^'*™'^ in terms of non-normal ordered condensates 

Uf-{q) = I ^ (vac|cc|vac) + -1^ (vacl^G^ac) - | ^ f In 4 + ^ V 
2 24 g"' vr 8 \ / 

(2.41) 

where we also have used the projection of the Lorentz structure for the vacuum 
case: (vac|c7Ai -D^iC|vac) = —^g\^ (vac|cc|vac); cf. Eq. flEllI) for the corresponding 
in-medium projection. 

Now we apply the heavy-quark mass (HQM) expansion, which is an expansion 
of a heavy quark condensate in inverse powers of the heavy quark mass, and for the 



scalar charmed condensate in vacuum it is given by 



36 
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53]: 



(vaclcclvac) = -ml f In At + 1^ ^ (vacl— CVac) + O (mZ^) . (2.42) 

4 V m2 J 12mc vr \ c / y 

The idea behind the HQM expansion is, that the interaction of the heavy quarks with 
the QCD vacuum mainly happens via gluon interactions, because the quark itself is 
too heavy to couple directly to a condensate. Hence, the heavy quark condensates 
are expressed in terms of gluon condensates. Obviously, by inserting (I2.42p into 
(I2.4ip we recognize that both the additional gluon condensate and the logarithmic 
term which appear in (12.411) are exactly cancelled by the HQM expansion. This 
statement of cancellation can also be found e.g. in 22|. 

Let us now consider the nucleon part: Hc^''^. The charmed condensates in a 



nucleon are negligibly small 



35| . since they contribute in nucleon only via virtual 
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effects; in this respect, to generate an easily interpreted image we recall the Tamm- 



Dancoff approach 



381-42 



-42] where the nucleon consists of a valence quark core of up 
and down quarks accompanied by a cloud of virtual mesons which accounts for the 
virtual sea quarks (e.g. charm quarks) and gluons. Thus, we can neglect the in- 
medium charmed condensates. Finally we note, that it is almost certain for us that 
the described cancellation of charmed condensates in vacuum via HQM expansion 
can be generalized to the in-medium charmed condensates, but a detailed proof 
needs special care and would spoil the intention of our paper. 

In summary of this Section, the complete OPE up to operators of mass dimension 
5 and up to the first non-trivial order in the coupling constant in terms of normal- 
ordered operators is given by Eq. (12.1 7p . where (I2.18P is given by Eqs. (I2.25P - (I2.3ip . 
and fl219|) is given by Eq. f l238|) : the term (^20^ is given by Eq. f l239|) but will be 
neglected in our further investigation. Finally, the term (I2.2ip does not contribute 
at large |g^|. 

III. OPE IN TERMS OF NON-NORMAL-ORDERED OPERATORS 
A. Absorption of IR divergences 

In order to perform a consistent separation of scales, all the infrared divergences 
have to be absorbed into the condensates, which means that the coefficient functions 
are only determined by the short distance behavior, while the non-perturbative 



effects are encoded in the condensates. In 



50-52] it has been shown that the Wilson 



coefficients are polynomial functions of the mass only when they are calculated 
in a minimal subtraction scheme. In contrast, normal-ordering is not a minimal 
subtraction scheme. Especially, by introducing normal-ordered condensates the IR- 
divergent terms appear explicitly in the Wilson coefficients. For a small quark mass 
these IR-divergent terms are remnants of the long-distance part of the correlator 
and have to be absorbed into the condensates since they basically contain the long- 
distance part of the OPE. Therefore, we have to express the normal-ordered tree-level 
condensates (: O :)''°^ in terms of non-normal-ordered condensates on one-loop level 
{0)^^\ That means, the OPE in Eq. (12. 2 p is reexpressed in terms of non-normal- 
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ordered operators: 

n(g) = ca+Y,c^{n\o,\n)^^K (3.1) 

i 

In the expansion (13. ip a proper separation of short and long distance contributions 
has been performed. Especially, the Wilson coefficients will be finite in the limit of a 
vanishing light-quark mass — )■ 0, that means IR-stable. The Wilson coefficients 
in (13. ip have no bar, which indicates that they are Wilson coefficients of non-normal- 
ordered operators on one-loop level. In consistency with (12.21) . we have to take into 
account all operators up to mass dimension 5. 



B. Operator mixing 



In order to express the normal-ordered condensates by physical condensates, we 
note Wick's theorem for an equal-time operator product of two quark fields. In 
doing so we partly follow |53| to get 



Tg(x)0 q{y) = : q{x)o\^^ q{y) : -i : Trc,D [o S{y,x)) : . 

(3.2) 



Here, O 



1^, 



denotes a function of covariant derivatives in respect to y, gluon fields, 
and Dirac structures. Now we set x = in the expression (13. 2p and insert a Fourier 
representation of the quark propagator: 



Tg(0)C 



3, 



q(y) = : m O 
(2 7r)^ 



q{y) : 
: Trcn ( O 



(3.3) 



Subsequently we take the limit y — )■ 0, i.e. no time-ordering anymore, and obtain 
the important relation 



{n\qO 



q\9) = {Q\ : qO 



q:\Q) 



{n\ : Trc,D O -ip, - S{p) ) : \n) ,(3.4) 



where we have also taken the expectation value. So far, the relation (13. 4p has been 



applied to the in-medium case in [33|, 



54 



. Here, O 



-ip^ - lA^ 



denotes 



21 



the Fourier transformed operator function, where Afj_ is defined by Eq. f IBlip . The 



derivatives are now contained in the gluon fields. Therefore, the ordering among the 
Lorentz indices is important as the field operators act on everything to the right. 
We define (13. 4p as the relation between physical condensates and normal-ordered 
condensates. Renormalizing the physical condensates by absorbing the divergent 
terms of 11^] on the r.h.s. of Eq. (13. 4 p cancels out the infrared divergences and yields 
infrared stable Wilson coefficients. We note that the extension of operator mixing 
for the four-quark condensates can straightforward be obtained by very similar steps 
like (1321) - (M . 

Eq. (13. 4p is valid to any loop-order of the corresponding Feynman diagrams. 
Here, for our purposes it will be sufficient to calculate the matrix elements only up 
to one-loop order, denoted by the index Up to order O (a^) the following set of 
relations has been obtained in MS scheme: 

{Q\ : ^ : (3.5) 



12m g IT 



(3.6) 



{^\q^^i^,qW^^ = {n\:q^,tt^q:\ny'^ + -^m^gg,J\n^ + l 



48 ^ ' TT ' ^ 



{n\qtt^tt,q\n)^'^ = {n\:qtt^tt,q:\n)^''^ + ^g^J\n^ + l] 

-i(«-^^)(-^n)<"|^vf¥4)^|")™ <3.8, 



In Appendix [D] relation (13. 8p is proven in some detail. Furthermore, we notice the 
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relations 

G'2|^])(i) = (^n\:^G': |^])(°) , (3.9) 

TT TT 

because these operators are already of order 0{as). The self-consistency of the 
relations (13. 7p and (13. 8 p with the relations (13 .Sp and (I3.6P can easily be verified by a 
contraction with the metric tensor and with the aid of the equation of motion ( ]A3p 
and relation ( lASp . respectively. 




The relations for the scalar operators (13. 5p and (13.61) are identical 



;o the relations 



for the vacuum condensates and were eiven, for instance, in |25|. l28l. 144 1: in vacuum 

n 

higher orders of (13.50 were obtained in [29\ . For the other condensates which vanish 
in vacuum we obtain of course {q^jf^q)''^^ = (: ^7/^1? ■)^^\ {qi ^q)*^^"^ = {:qi ^q '■)^^\ 
and {q'^->laG^uqY^'' = {■ qiblcxG q ■.)^^\ 

From now on we drop the label i.e. non- normal-ordered condensates are 
always on one-loop-level throughout the work. We note, that a scale /i naturally 
appears when the non-normal-ordered condensates are introduced, and they would 
even not be well-defined without /i. The /i scale is an IR cutoff which marks the 
region between perturbative and non-perturbative physics. Hence, any sensitivity 
to the light-quark mass must be accounted for by a redefinition of the condensates. 
This resembles the case of parton distribution functions (PDFs) in deep inelastic 
scattering. The next-to-leading-order calculations are IR divergent which can be 
accounted for by a redefined PDF. As a consequence the PDFs become scale depen- 
dent which gives rise to the famous DGLAP (Dokshitzer, Gribov, Lipatov, Altarelli, 
Parisi) equations 59l-l6l|. We underline that a careful application of these operator 
mixing relations (13. 5 p - (I3.10p is mandatory for a consistent treatment of the OPE if 
one internal quark mass of the currents is much smaller than the confinement scale 
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C. Insertion of operator mixing 



By inserting Eqs. (13. 5p - (13.101) into Eq. (I2.38P we obtain 11^^^ (g) in terms of 
non- normal-ordered condensates: 



uf\q) = {n\dd\n) 



rrir 



q2 _ ^2 



-{n 

-{Q 
-{Q 



q 



A 



m 



2\3 



m 



2^2 



TV 



1 nic 



1 



TT 



[vG] 



12 md g2 

2 ^ 



24 (g2 
1 



2^2 



9 (g2 



2^2 



g2 -4 



24 g2 



In 



-77 K-S-ll) 



For the mixed quark-gluon condensate (next to last term in Eq. (I2.38P ) we have used 
relation ( lASp . while for the mixed quark-gluon condensate which contains 75 (last 
term in Eq. (12.380 ) we have used the relation flE22p which can be shown to be valid 
for an arbitrary four-vector. Special care is needed when applying the equation of 
motion flA3p to some condensates in (I2.38P and neglecting terms oc m^, because the 
operator mixing may introduce terms which cancel out factors of m^. For instance, 
this is the case for (: d'yxl5 ^ :); contraction with the metric tensor and application 
of equation of motion flA3p results in a term oc md{'- dd :), which, however, remains 
finite and non-zero after operator mixing. Thus, we cannot neglect such terms. 

Here, we emphasize that the both IR-divergent terms oc m^^ and oc Inm^ in 
Eq. (13. lip will cancel the corresponding IR-divergent terms in (12.3 ip . This cancel- 
lation will be the topic of the next Section. 



IV. OPE IN THE LIMIT OF A VANISHING LIGHT-QUARK MASS 

A sensitivity to light-quark mass signals a sensitivity to very small quark and 
gluon momenta. Physically such momenta are cut off by confinement, i.e. at a scale 
of Aqcd — 245 MeV, and not by the light current quark mass at a few MeV scale. 
Hence, any sensitivity to the quark mass must be accounted for by a redefinition of 
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the condensates. Moreover, the hmit — )■ is important, since it demonstrates 
the cancellation of all IR-divergences in fl2.3ip and f l3.1ip . which is a fundamental 
feature for the consistency of the whole approach. 

In order to calculate the complete OPE for the D-meson in medium with infrared 
stable Wilson coefficients in terms of non-normal-ordered condensates we have to add 
the perturbative part (12.251) . the expressions (I2.3ip by virtue of Eqs. (13. 9 p and (I3.10p 
and the terms in Eq. (I3.1ip . We emphasize the cancellation of the IR-divergent term 
oc In(m^) in the Wilson coefficient of the medium- specific gluon condensate and the 
IR-divergent term oc m^^ in the Wilson coefficient of the gluon condensate. This 
important result of cancellation of IR-divergences due operator mixing (13. 4p has 



been demonstrated at the first time in 



54| for in- medium D- mesons at rest. In a 



more comprehensive analysis this result has been extended to the case of in-medium 
D- mesons with finite three-momenta in 55 1. 



Now we will take the limit — )■ which shows explicitly that the complete 
OPE in terms of non-normal-ordered condensates is infrared stable. The complete 
projections with respect to Dirac- and color- indices of the condensates onto invariant 
structures are given in Appendix [El By means of these projections, we obtain up to 
mass dimension 5: 



n(g) = ^p^'-*(g) + (^]|(Jrf|^])- 
-(^]|^G'2|^]) ^ ^ 



{VL\dgagd\9)\ "^""^^ 



2 (g2 



2^3 



X - 




{vqf 



54 (g2 



m 



2^2 



9 (g2 



m 



2^2 



In 



9 (g2 



m 



2^2 



In 



.(^Q\dv^l^d\Q){q^ 



{vqf 



3 iq' 



m 



2^2 



-{n\dv{v'tYd\n)^ (q^-2 



Xvqf 



+{n\dvgs(Tgd\n) 




(g2 



m 



2^3 



(g2 

{n\dgsagd\n) 



m 



2^2 



{Q\dvd\n) 



vq 



y2 q2 



mt 



{vqf 



rrir. 



(g2 — m^ 



|3 ■ 



(4.1) 
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The angled brackets in the last term indicate that this combination vanishes in 
vacuum, thus denotes a medium-specific part: applying vacuum projections to this 
medium specific term makes it zero in the vacuum limit. Numerical values of the 
condensates in f l4.ip are given in 33j. The IR-stable OPE in Eq. f l4.ip is the main 
result of our investigation. It is valid for a meson whith the four-momentum = 
{qo,q), while the surrounding medium has a four- velocity = {l,v) in respect to 
a given frame. 

Needless to say, in the limit — )■ the OPE simplifies considerably. But 
we emphasize that the limit ma — )■ is not performed because such terms are small 
anyway or in order to simplify the OPE, instead this limit is necessary from physical 
reasons: the dependence of the correlator from the light-quark mass m^; signals a 
sensitivity to very small momenta, which are cut off by confinement, that means 
at a scale of Aqcd, which is much larger than the light-quark mass. Thus only the 
OPE in the limit of a vanishing light-quark mass is meaningful from the physical 
point of view. 



V. OPE FOR A MESON AT REST WITH RESPECT TO THE MEDIUM 

In the previous Section we have presented the OPE for the general case g'^ = 
{qo,q)- However, for many investigations concerning in-medium properties of D- 
mesons it will be sufficient to consider the D-mesons at rest with respect to the 
medium they are embedded in. Therefore, we will also consider the special case = 
(go, 0); furthermore we choose a frame comoving with the medium v'^ = (1, 0). We 
present the OPE separated into an even part and an odd part: n(go) = n'^™"(gQ) -|- 
go n°*^'^(gQ). Furthermore, we perform an analytical continuation go = iu. Then, 
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the even part of the IR-stable OPE for — > can be written in the form: 

— 771 — 1 771 iiP' 



12 TT u;^ + 



2 " 2 

-2{Q\d^it)Qd\^) ^ ^ ^ - 4 {^\dt)ld\^) - {n\-dgsagd\n) 



{u^ + ml) 



mcUJ 



(5.1) 



Note, the angled brackets in the last term denote that this combination vanishes in 
vacuum, thus denotes a medium-specific part which is absent in vacuum. The odd 
part of the IR-stable OPE for — > is given by: 



TT^^iu) = {n\Sd\Q)— - + A{n\S^ld\Q) 



u:'' + ml ' ' ° ' ' (a;2 + mlf 

-{n\d^g,agd\n) ^ (5.2) 

[u'^ -\- m^) 

The OPE is an asymptotic scries. As such it must be truncated at a certain mass 
dimension and, therefore, takes into account a finite number of operators. One way 
to deal with asymptotic series is to perform a Borel transformation which suppresses 
the effect of higher mass dimensional operators. The Borel transformation of a 
function /(Q^) is defined by 

B [f (A'P)] = hm hm ( '-LX f iQ') , (5-3) 

where the parameter M is the so-called Borel mass. Applying the Borel transfor- 
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mation to the even part (15. ip we obtain 

oo 

B [n"™"(w2)j = - J ds e-'/^' ImnP^'-*(s) 



+e" 



-ml/M^ 



mc{^\dd\^l) + - 



1 / mi 



TTlr 



1 



, , ^]{Q\dgsagd\Q) + — (Q\ — G^\Q) 



+ 



T_ 1 /i^ mg _ 2 ...^ 
18 ^3 ""ll^ 3^^/ VM2 



2 

3 M2 



/ {vGf G 



-2 (^-l) (^M^^^orfl^^) 



+4 



2M4 ~ 



(5.4) 



while for the odd part (15. 2p we obtain 



2M4 M2 



{Q\d^ 'dld\Q) 



M2 



(5.5) 



This Borel transformed OPE has been presented in 33|, |3J| for D- mesons in-medium, 
however no further details about how to arrive at this correct OPE have been pre- 
sented; so far a detailed presentation can only be found in js^l and in the more 



extensive study where the a ppr oac. 



divergences for the vacuum case 25 
case. 



1 of operator mixing and cancellation of IR- 
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29| has been generalized to the in-medium 



VI. COMPARISON WITH THE LITERATURE 



Even the OPE for D-mesons in vacuum is not as trivial as the common belief 
might be. This can be illustrated, for instance, by a brief review of the history 
about the different results obtained for the Wilson coefficient of the condensate 



(vac|(7s (io" ^ (i|vac). The first attempt to calculate this coefficient was done in 62l |. 
where a wrong factor —1/4 has been presented. At the first time, the correct result 
—1/2 has been presented in 20 1, where the needed correction of |62| was explicitly 



mentioned. Later, in 



63| a wrong factor +1/4 has been given. The correct result in 



28 



Re f. [20 1 has later been confirmed in 28, 64|. But this was not the end of the story. 
In 65| a wrong factor has been given again, which has later been corrected in |66| 

65|. 
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by the same author, exphcitly mentioning the needed corrections in Refs. 
This brief survey of history shows, that the calculation of the OPE for D-mesons 
needs special care, and even more for the in-medium case. 

So far, in- medium QSR for D-mesons were given in 30N32| and in our own inves- 



tigations 



l55l. The references 



3l| and 32| are from the same author(s) 



and the given results agree with each other. Furthermore, in a recent study |64| the 
D-mesons in vacuum were considered, but since the Wilson coefficients of scalar con- 
densates in vacuum and in-medium are the same we can compare the results. Thus, 
at the moment being a comparison of our results is meaningful with the findings of 
Refs. laO, 31, 64|. 



A. Scalar part of the OPE for D-mesons in medium 



In Table [T] we compare our results for scalar condensates in-medium with the 
corresponding results of Refs. 



30 



31 



64j . Here, we list the Borel transformed Wilson 



coefficients c^/e 



-m2/M2 



Obviously, there are serious disagreements of the findings 



in [30|, |64] to our results, while the results of Ref. [3J| do agree. It seemed to us, 
that the wrong Wilson coefficient for the scalar gluon condensate in 30|, l64] is either 
caused by a missing sign in an intermediate step somewhere in their calculation or 



because of an incorrect operator mixing in Refs. 30|, |64| ]. 
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11/1 1 rn'i) 
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TABLE I: Comparison of scalar condensates of the OPE given by Eq. (|5.4p with Refs. 
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ej. In 



3C 



301 ] there was no scalar mixed condensate, because only operators up to mass 



dimension 4 were taken into account. 
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B. Tensor part of the OPE for D-mesons in medium 



In Table HTl we compare the tensor part of the in-medium OPE with the litera- 
ture. As before, we list the Borel transformed Wilson coefficients co/e"™^/^^. We 
have found a serious disagreement between our results and the findings of Ref. [30 1. 



Un: 
in 



brtunately, it is impossible to compare our results directly with the ones given 



31 



32|, as it is somehow hidden in their work. We could not understand in 



detail how the authors of [311, l32| have treated the IR singularities. However, for the 
unprojected OPE and before introducing physical condensates such a comparison 
is possible. That means we compare our result (12.381) with 3l(|, although this is 
less meaningful, because many calculations still have to be done from this point on. 
We have found an agreement of fl2.38p with jsi] to a large extent, except a slight 
disagreement of a factor | in the Wilson coefficient of the condensate (: d'jxaQd :), 
which emerges from (: dl^ i^yd :) and (: d'y^^aQfiud :) by means of the equation of 
motion (1A3P : for further details see Ssl . 





{Vt\dU^od\VL) 




Eq. ([OD 
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14 I ml ,1 ml (. ml\, ( ml ^ 
3 3 eJP^2W^^\^ Jl^ ) '■^^ \Atv ^i^ J 

+^ml/M^l 2 7ij ln^+ / dt^-f\ 
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TABLE II: Comparison of tensorial condensates of OPE given by Eq. (|5.4p with Ref. 



Finally, we note that for the odd part of the in-medium OPE a comparison of our 
result with [30]] is also not possible, because it is not clear whether the odd part 

jart in vacuum, hence a 



has been considered in l30l or not; note there is no ode 



32| the OPE is given in 



comparison with [64] is impossible. Furthermore, in [31 
terms of normal-ordered condensates which implies IR-divergent Wilson coefficients 
and makes a comparison difficult. 

In summary, we come to the conclusion that the OPE for in-medium D-mesons 
has not been determined correctly, except by the presentations given in 



33 



34| 



where, however, no further details of the involved calculations have been given; so 



30 



far, such details were only presented in 



54| and in the more detailed analysis 



55| 



VII. SUMMARY 

We have determined the in-medium OPE, at first non-trivial order in powers of 
the QCD coupling constant and including all operators up to mass dimension 5, for 
heavy-light-quark pseudoscalar mesons: D-mesons and B-mesons. We have outlined 
that it is sufficient to consider the OPE just for mesons, since all other OPE's 
of heavy-light-quark pseudoscalar mesons can easily be deduced from that result. 

So far, in-medium QSR for D-mesons are fairly rare 
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54i5;l and the OPE 
32| is incorrect. 



in these works differ significantly. Especially, the applied OPE in 130 
From our view, the reason for this fact is that the derivation of the correct OPE 
turns out to be an ambitious assignment of a task. Accordingly, it is timely to 
present the derivation of the OPE for D-mesons in matter in some detail. 

First, we have determined the OPE with the aid of Wick's theorem, leading to an 
OPE in terms of normal-ordered condensates, see Eq. (12.21) . The complete OPE up 
to operators of mass dimension 5 in terms of normal-ordered condensates is given 
by Eq. (l2Tn) and Eq. fTM the term by Eqs. ^22^ and (|2:26|) . the term Y\!§, 
by Eq. fl2.3ip and the term 11^^ by Eq. fl2.38l) . We also have determined Ili^'' given 
by Eq. fl2.39p which contains charmed condensates, but we have argued why these 
charmed condensates can be neglected. The Wilson coefficients of OPE (12. 2p are 
IR-divergent, that means they are infinite in the limit of a vanishing light-quark 
mass, md — )■ 0. 

It has been described in detail that a consistent treatment of the OPE is ob- 
tained in terms of non-normal-ordered condensates, see Eq. (13. ip . The relation 
between normal-ordered condensates and non-normal-ordered condensates is given 
by Eq. (13.41) to any loop-order. We have determined explicitly the relations among 
these condensates to one-loop-order, that is nothing else but the operator mixing 
under finite renormalization, see Eqs. (13. 5 p - (I3.10p . By means of these relations we 
have obtained the OPE with IR-stable Wilson coefficients which allows the limit of 
a vanishing light-quark mass. 

The result of an OPE in the limit — )■ is given by Eq. (14. ip . which is valid for 
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moving D-mesons with respect to the surrounding medium, first obtained in 55 1. 
The OPE in Eq. (14. ip is the main resuh of our investigation. The important case of 
D-mesons at rest with respect to the medium simphfies the OPE considerably and 
is given by Eqs. f lS.ip and fl5.2p . and has first been obtained in 54 1. 

The OPE for a D-meson at rest has been compared with other OPE's used so far 
in the hterature. Remarkably, we have found that some theoretical investigations 
have used seriously incorrect expressions for the OPE of D-mesons in matter. The 
aim of our investigation is, therefore, to present a more detailed analysis about how 
to obtain the correct OPE of D-mesons in medium. We hope that our investigation 
will support the correctness of prospective theoretical studies. 
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Appendix A: Notation and conventions 

In this Appendix we briefly define the basic quantities and conventions used 



throughout this work. Our expressions are obtained within the framework of [67 1. 
Contravariant four- vectors = {qo,q), covariant four- vectors = {qo,—Q) in 
Minkowski space. The Lagrangian density of QCD reads 
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where M = diag (m„, m^, m^, mc, m;,, mf) is the matrix of current quark mass, the 
quark field flavors are \E' = {u, d, s, c, b, t)'^, and the covariant derivative 



(A2) 



with Qs = ^^nag being the strong couphng constant, and X°-^ is the unit matrix 
in color space. The notation means not a vector but denotes that the partial 
derivative acts on the field to the right, vice versa acts on the field to the left. 
The gluon fields are J^f^ = t'^f^ and the gluon field strength tensor Q^^ = G^^,t^j^; 
recall the use of Einstein convention. Here, = | Af^ are the generators of SU (Nc) 
where A^^^ = 3 is the number of colors and A^^ are the Gell-Mann matrices. The 
Lorentz indices are denoted by greek indices /i, = 0,1,2,3, the color indices are 
denoted by small letters a,b = 1,2,3, and the Gell-Mann index is denoted by a 
capital letter A = 1,...,8. Finally, we note the equation of motion for quark and 
gluon fields: 



—I niq q 



/ 



(A3) 
(A4) 



We also note the needed relation 



m 



(A5) 



where 



gf"" Uy. i-urtnermore, we use a^^ = ^ [7^ , = i (^7^7^ - g^,^}, 
where the metric tensor is g^y = diag (+1, — 1, — 1, — 1), and for the Levi-Civita 



Furthermore, we use cr^,^ = f [7/, , 7v] 



symbol we use the convention e 



afiucr 



afii/ 



Gg''^ with ei234 = +1- 



Appendix B: Background field method in Fock-Schwinger gauge 



The OPE is calculated using the background field method 22[ by applying Wick's 
theorem to the current-current correlator and employing the Fock-Schwinger gauge 



{x^-x^^)A,ix) 



0, 



(Bl) 



first introduced in Quantum 



later extended to QCD 



22 



01, 



Electrodynamics [68|, l69| for the photon fields, and 
70| for the gluon fields. Usually one choses the frame 
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Xq = 0. Very detailed aspects about translational invariance and the employed 
gauge can especially be found in [6^]. One decisive advantage of Fock-Schwinger 
gauge is, that the gluon fields can be expressed in terms of the gluon field strength 
tensor: 

1 





This relation allows for a considerable simplification of the calculations. Further- 
more, the Fock-Schwinger gauge allows to express partial derivatives of quark and 
gluon fields by covariant derivatives. Accordingly, a covariant expansion for gluon 
and quark fields in Fock-Schwinger gauge is given by 

oo y 

fc=0 ^ ' 
oo ^ 

^(^) = E fc!^"' • • • (^-1 • • • ^°'^'^).=o ' ^^^^ 

k=0 ' ^ 

oo _ 

(^^r^„^...^^^j^^^ . (B5) 



fe=0 

The propagator of a quark in a weak gluonic background field in coordinate space 
is given by Eq. (E^S]). 

We will now give a closed expression for the quark propagator in momentum 
space which can be calculated directly from Eq. f l2.23p by introducing the Fourier 
transforms of the quark propagator and the gluon field: 

S{p) = j d'^xe'P'' S{x,0) , S{p) = j d^xe-'P''S{0,x), (B6) 

A^ip) = J d^x e'P^ A^{x) . (B7) 

For the Fourier transform of gluon fields flB7l) we obtain 

(B8) 



where the partial derivatives are with respect to four-momentum. By repeated 
partial integration and careful consideration of partial derivatives acting on Dirac's 
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delta distribution, one can confirm the following form for the quark propagator 
Eq. fl2.23p in momentum space 



S{p) = (B9) 
= (_i)'=5(o)(p) (^^ij 5(0) (p) X ... X (7^) , (BIO) 



k 



where A is a derivative operator which naturally emerges during the Fourier trans- 
form and is defined as 

00 

k=0 

= (^"1 • • • d^d-^ ...d-^. (B12) 

Due to the Fock-Schwinger gauge, the quark propagator (12.231) does not obey trans- 
lational invariance and, hence, is not a function of the coordinate differences and 



cannot be written as S{x, y) = S{x — y) 67|]. However, after performing the Fourier 



transformation, the difference between S{x,0) and S{0,y) is merely manifested by 
the operators acting on the terms to the right or to the left of them and the 
equality can be read off. 

For the OPE up to mass dimension 5, only the first three terms in the sum of 
Eq. (1B9I) became relevant, and are given by {n is either c or d) 

si^\p) = 4^^'^^ (B13) 

S^r^\p) = ^^QtMS'^\p)rs'^\p)^^S^\p), (B14) 

S^\p) = I)' ^;^(0) ^:a(0) Tr\p) , (B15) 
where we have defined the tensor 

Tr\p) = sl^\p)YS^:\p)Ysl^\p)i''s^:\p)^' 

ip) r (p) 1- (p) Y (p) 7^ 

(p) 7^ S^:^ ip) r ip) 7' ip) Y ■ (B16) 
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Appendix C: The terms ([27281) " (12130]) 



The evaluation of the gluonic expressions fl2.28p - fl2.30p is involved and leads 
to intricate expressions. We therefore would like to give some basic steps about 
how to arrive at the final expression given by Eq. (I2.3ip . including the needed 
approximations. 

First, by inserting the expressions flB13P - flBlSP into f l2.28p - f l2.30p we obtain 



X (75 Si'^ ip) Si'^ ip) Y Si'^ (p)75 (p - q) Y Sf^ {p - q) Y Sf ip - q) 



(CI) 



2 



X / ^Tro (75 (p + q) 75 Tr\p)) , (C3) 

where the trace over color-indices has been performed, i.e. TrD denotes the trace 
over Dirac-indices only, which are evaluated by means of the algebraic software HIP 
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7l| . Finally, we are left with the following four basic integrals: 
d^p 1 



4-l,fc-l,n.+fc-2(g^ "^d, ml) , (C4) 



1 r(n + A;-2) ^ ,222 



(47r)2 T{n)T{k) 

d'^p p^ 



(2 7r)4[(p_g)2 + ^2]n [p2 + ^2]* 

In,k-i,n+k-2{q^ , m'i, ml) , (C5) 



T{n + k-2) / 2 2 2 



(4 7r)2 r(n)r(A;) 



(2 7r)4[(p_g)2 + ^2]n [p2 + ^: 



5f^j. 1 r{n + k-3) 



(4 7r)2 2 T{n)T{k) 



g^g^ r(n + A;-2) 3 2 2^ 

(4^ Tin)Tik) ^r.+i,k-i,n+k-2{q , m„ mj , (C6) 
d'^p PplPvPk 



(2 7r)4[(p_g)2 + ^2]n [p2 + ^2] 



Qi^u qn + 9^ + 5'^^. i r(n + A;-3) / 2 2 2n 

a^^q^qn T{n + k-2) 2 2 2 



(4 7r)2 r(n)r(A;) 



In+2,k-i,n+k-2{q ,ma,m^) , (C7) 



where r(x) = j dtt^ ^ e * is the gamma function, and the master integral reads 


1 

k^^^i^^, ml ml) = / da —- ' ' • (C8) 

J [a[ \ — a) q'^ + a + [l — a) 

Needless to say that these integrals lead to rather cumbersome expressions for finite 
mass rrid and m^. On the other side it is not possible to put them to zero since the 
integrals are IR-divergent, that means they contain terms like m^^ and Inm^i. In 
order to extract the corresponding terms cx m^^ one has to multiply the integral with 
the lowest power n of that gives a finite result for the mass going to zero. What 
remains is the corresponding coefficient for the m^" term. All the other coefficients 
can be obtained by taking suitable derivatives of m^Iijk with respect to and 
setting afterward equal to zero. The terms oc Inm^ can be obtained by taking 
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the first derivative of lijk and extracting tfie term oc m^^ from that expression in 
the same way. In this respect we also note the useful derivative relations 

d 



d 
drnl 



Iijk{q^,ml,ml) = -klij+i^k+i{q^,ml,ml) , (CIO) 



d 

—Iijk{q^,ml,ml) = -/c/i+ij+i,fc+i(g^ m^, m^) . (Cll) 

Sometimes, some additional manipulations have to be made in order to obtain mean- 
ingful expressions. Using 

1 1 + iz 

arctanz = — In (C12) 

2i 1 — iz 

one gets 



2 2 \ -I 

q — \ 1 m 



arctan 1 = in ^ , (C13) 

\ q + J 1% 

which is a source of mass logarithms. Another source of terms oc Inm^ arises from 
a slightly different expression, namely 

arctan [ i ^ q^ + ml-ml \ m^o ^^^^^^^^^ ^ ^^14) 

y y 2g^m^ + + 2q^m1 + — 2m^m^ + j 

which is not well defined. Expanding the fraction in and keeping only the lowest 
power, which is the dominant contribution for rud — ?■ 0, leads to 

q' + K-K « 



2q^ml + + 2q^ml + - 2mlmj + mj J 2i (g^ + ' 

(C15) 

This term arises in the tensor part of the term (ICSP and causes an infrared divergent 
Wilson coefficient. It prohibits us from taking the limit rrid — > at this stage. 
These formulas and the technique described here in detail allow us to evaluate the 
integrals (ICip - (lC3p . We underline that the expressions (ICip - (ICSp are symmetric 
in exchanging charm and down quarks. However, it is obvious that the needed series 
expansions in respect to the small quark mass destroy this symmetry. We finally 
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have obtained the following results: 



nf:&) = +^(^^|:fG^:|^^)^^, (C17) 

^ _ . 2 . ( I 1 1 1 _ 1 rnl \ 

(C18) 

where = G^J^^- We recognize that the term oc m^^ prohibits us from regarding 
one quark to be massless in evaluating the scalar terms. For the tensor part we obtain 
for rrid ~ the following results: 

ntensorl?) = (^M — — i ^ ) ■ \^) [ Q - 4 



f 1ml f ml \ 1 1 1 , 

n^^^r^ /Ol "^^M' (, Avqr\f 1 1 1 

ntensorl?) = (^h — 5 T ' 1^/ U " ^ 



7r\yf2 4y \ "^^yV -'-89^ 9^ ~ ""^c 



7r\yf2 4/ \ J \ 18g2 g2 _ ^2 

/ 1 2m2\ , / 



Qq"^ {q"^ — \Qq^ 9q^ J \ q^ — ml 

1 / mi 1 \ , /m' 



+ ^ In ^ 



9g2 \{^q^ — mlY q^ — mlJ \Tnl 
-— ( + ^ l In ( 1 1 . (C21) 

Here, we observe the occurrence of a term cx Inm^, which again prohibits us from 
taking the limit — )■ 0. The divergent terms cancel, however, after the introduction 
of physical condensates. According to Eq. f l2.27p . we collect all the terms fIClGp - 
( 1C2ip and obtain the expression given by Eq. fl2.3ip . 
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Appendix D: Proof of relation (13. 8|) 



The operator in relation (13. 8p reads O 



From the operator 



mixing f l3.4p we obtain in one-loop approximation 

(M 



(fi|:Trc,D {p, + Aj{p, + A,]S{p) :|fi)W.(Dl) 



In consistence with the whole approach, we determine the integral up to order 0{g'^), 
and the expression in the last line can be separated into three terms as follows: 

d^p 



T3 



—I 



(2vr)4 
d'^p 

d'^p 
(2^ 



(D2) 
(D3) 
(D4) 



while all other terms vanish. Let us consider explicitly Ti. Using the expressions 
of the fermion propagator given by Eqs. (]B13P - (IBlSp . and the derivative operator 
Q^y which acts on the quark propagators only, we obtain 



i(o) _ !^ 



-1^) (^^|:GiGf,:|r])(°)Trc 



2 2 



2e 



d^p 1 
(27r)^ (p2 _ ^2) 



gTro 



(D5) 



where the integrals are divergent and have to be evaluated in D = 4 — 2 e dimensions, 
i.e. by means of dimensional regularization; u is the renormalization scale of finite 

nn 

renormalization governed by the renormalization group equation p^, 119[ . 

Now we insert the projection (IE12I) and perform the Dirac traces; note the color 
trace is TrpA^ A^ = 2 . After a lengthly but straigtforward calculation we finally 
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obtain 



2e 







2 


) 96 






/ 


(27r)« 



X/i 



(24 ml + 8mqp'^ g^y - 32 



2 ^2 

— 1 : 



4 



— « 



2e 



d^p 1 



(24 + 8 mg gf^,, - 32 p^p^) 



2t 





y 8 


2 


y 24 




d^p 


/ 





v2 4 



(0) 



^ -3 ( -Snigip'^ - ml)g^y + IQmqP^Vf.v^ + 16 



-32mq {pv){v^pu + VuP^,) + 32 mq{pvfg^u 



mq v^Vy 



(D6) 



where we have kept some ratios in order to show more exphcitly their source from 
the projection. The two needed integrals can be evaluated by standard techniques 
and in MS scheme 



47( 1 they are given by: 



1 1 



d'^p 1 
(2^ {p2-mlY~'2^gj4^ 

1 

~4 



2e / PlJ^P. 



1 i 



l-ln4 

2 m2 



(D7) 
(D8) 



(27r)^ (p2_^2)3 

The second integral is logarithmically divergent and has to be evaluated by means 
of dimensional regularization. Inserting both integrals into (lD6p . we arrive at 



^1 = +^mqg,,{n\:'^G^:\n) 
6Z n 



(0) 



24 



4 



(D9) 



Now we consider the term T2. Using = 1, and inserting (1B13P into (1D3P yields 

d^p p^ 



T2 = -3 ruq gf^iyijj'^^ 
41 



(27r)^ p2 — rriq 



(DIO) 



where we have used that the only possible Lorentz structure is the metric tensor, 
and we have used Trc [S'^^~\ = 3, TrD [p] = 0, TrD [rriq] = Dniq and (7^' = D, where 
D = 4— 2 e. The nominator of the integral is treated by means of = {p^ —m?^+m^^, 
and then using the common law j ^^^^ = by definition, e.g. 73|, we conclude 
the superficial degree of divergence of this integral is 2. Standard evaluation in MS 



scheme 



47| yields: 



Thus, we obtain 



2e 



{27r)^ p'^ — 



(4vr)2 



lni^ + 1 



To 



3 

16 TT' 



;9f,u In 



ml 



+ 1 



(Dll) 



(D12) 



The determination of T3 is a lenghtly calculation because of the more involved Dirac 
traces, so we present just the final result. An application of the same techniques as 
used by Ti yields for T3: 



5 

96 



16 mi 



+ 



mi 



216 36 



5 

54 



mi 



mqV^Vy : 




(D13) 



According to Eq. (ID II) we add the terms Ti, T2 and T3 and finally obtain: 



{\nqi^^ii,q\n)^^'^ = {n\ : qi^ : \n) 



(0) 



+ 



m„ / , M 



16 

m 



{n\ : — : 

TT 



(0) 



9 / . . ^M^i^ 



36 



9ti,y - 4- 



3 



TT 



3 

16^ 



( In — + 1 



: |fi)(°),(D14) 



which is just relation 



Appendix E: Projections 



Condensates are vacuum or in-medium expectation values of quantum field op- 
erators. They reflect basic properties of the ground state of QCD or of the medium. 
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Condensates are assumed to be color singlets, Lorentz invariants and invariants un- 
der parity transformations and time reversal, according to the corresponding sym- 
metries assumed for the medium. Therefore, one has to project out color, spinor 
and Lorentz indices from several structures that appear during our calculations. 
Expectation values which are not invariant under parity transformations and time 
reversal are supposed to be zero. In what follows we adopt the method described in 
72|. 



1. Projection of color and Dirac structure 

Up to mass dimension 5 we meet the following structures 

{mq'l^) , {qi'^,)%^\n) , {n\ (q^.'^.y q]\n) , (El) 

where the Dirac indices i, j = 0, 1, 2, 3. The last structure is already invariant under 
color rotations, thus one does not need to take care of color projections for this one. 
One can expand the remaining other condensates using an orthogonal basis. For 
Nc X Nc dimensional matrices such a set is given by the generators tj^^ of SU{Nc) 
supplemented by the unit matrix X in A^^^ dimensions; an appropriate scalar product 
is given by the trace operation {A,B) = Tt{AB). Using SU{3) generator identities 
the prescription for the expansion in terms of the generators of SU(A'^c) reads as 
follows 

Mab = ^CAtt^^CxXah , (E3) 

A 

CA = Trc(t^t-^) ^"^^^^^^ ^^"^ cx=^Trc(M). (E4) 

In order to retain color singlets the only non-vanishing expansion coefficient is Cx- 
All the other coefficients belong to expectation values which do not transform as 
scalars under color transformations. 

The projection of spinor indices proceeds in an analog way. A complete set is 
given by the elements Ok of the Clifford algebra, i.e. Ok G {X, 7^, cr^i^, Z757^, 75}, 
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satisfying Tid (OiO^) = iSf. The expansion reads 

{mO^,...qj\n) = 5^40j,, dk = ^{n\qOkO,,,„q\Q), (E5) 
k 

resulting in 

{n\q,o,,...qj\n) = ^(^{n\qO,,...q\mJ^ + ml,o,,...q\n)J^, 

+^{Q\qa^^Of,^,„q\Q)a!^^ - {Q\q^5'y^O^^,„q\Q)^5^!^i 

+ {n\q-f,0^,,„q\n)j5j^. (E6) 

The notation used is: 757^*^ = (757^)1^ and 75^4 = (75)1^. Together with color 
projection one gets 

{mO,,...q'\Q) = —( {Q\qO,,...q\mj^+ mi,0,,...<l\m 



+ ^{Q\qcx^^O^^,„q\Q)a^i - (yL\q-i^-i^0^y,„q\Vt)-i^-i^^ 

+ (fi|g750^....g|fi)75,,i). (E7) 



Here C^j,... stands for an arbitrary operator with Lorentz indices fiu ■ ■ ■ . In Ref. [72 
it has been stated that terms corresponding to the projection onto cr^j/, 757^, 75 do 
not appear due to parity and/or time reversal. This statement can be misunder- 
stood. It is true for 757^, 75, however not for a^^,, which does not contribute because 
there is no independent Lorentz structure that reflects the symmetry properties of 
a^u- In fact, it contributes if the gluon field strength tensor enters the operator 
product (i.e. O^^ = Q^y). 



2. Projection of Lorentz structure 

A striking difference of vacuum and in-medium projections appears when project- 
ing Lorentz indices. It is important to note that this projection is not an expansion 
in terms of a complete orthogonal set. It strongly depends on the structures avail- 
able to perform the projection. In vacuum, there are only two independent objects, 
the metric tensor g^^ and the total antisymmetric symbol e^yai3, being a pseudo- 
tensor under parity transformations. In medium the condensates also depend on 
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the medium's four-velocity which, therefore, is an additional structure for pro- 
jections. As a result, also pseudo-vectorial structures can be invariant under parity 
transformations. 

Now we give a list of the in-medium projections up to mass dimension 5. Terms 
that violate time reversal or parity invariance are omitted. Equations of motion 
enable us to rewrite the condensates in terms of canonical condensates 53 1. 



a. Condensates of mass dimension 3 

The in-medium projections of the Lorentz structure of the condensates of mass 
dimension 3 read: 

{n\qq\n) = {n\qq\n) , (E8) 
{n\q^,q\n) = {n\qvq\n)^. (E9) 

A condensate of the type (qcr^uq) occurs neither in vacuum nor in a medium since 
there is no possibility to create an antisymmetric structure in the Lorentz indices 
/i, z/. Condensates of the type {q'j^'j^iq), {q'Jnq) can not be projected onto structures 
that are invariant under parity transformations. 



b. Condensates of mass dimension 4 

The in- medium projections of the Lorentz structure of the condensates of mass 
dimension 4 read: 



{n\q'^^q\n) = -{n\qvq\n)i 



rriqVf, 



{n\q^,D,q\n) = -{n\qq\n)t^g,, 

i'I^{n\qq\n) + {n\qv^-^q\n) 



3 i^^" ^ y 



S 



AB 



96 



24 



(ElO) 



(Ell) 



(E12) 
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where we have defined 



a _ f of ^^'"^ ^/^^^ ^ 



3v 9 
V 



VaVu 

9/3 ij, — 5- 



(E13) 



Again, terms of the form (f?757^^^5'), {qI^^hQ) do not have a projection due to the 
requirement of parity invariance. The same holds true for {q<J ^iJl^ kQ) ■, which can be 
contracted with e^^^xv'^ giving an odd term with respect to parity. 



c. Condensates of mass dimension 5 



The in-medium projections of the Lorentz structure of the condensates of mass 
dimension 5 read: 

o 

{n\qD^D,q\n) = -{n\qq\n)^g^,+ {n\qg,agq\n)-g^. 



m 



1 



-^{n\qq\n) - -{n\qg,agq\n) + {n\q 



_{vDf 



4 8 



q\Q) 



(E14) 



{n\q-f^t,t^q\n) 

1 



^{fl\qv{vDyq\fl) ( ^ ''^7''" ~ \ {'"l^9ua + v„g^a + Vag,j.u)^ 



1 



,2 



{fl\qvgsagq\fl) 



2 Vi^V^Va 1 



, "^9 /oi- IO\ f^t^^^'^^ 
6V \ V 



(E15) 



(f^kTsTa^Mi^gl^^) = --^{VL\qvgs(Tgq\Q)eanuav'' , 



(E16) 



(fi|g5(3(T„/3^/,^g|l^) = {VL\qgs(7gq\Vt)— {g^^gf^^ - gau9f3^) 



+ 



1 2m (vD) 2 

-{n\qg,agq\n) -i^{Q\qv^^q\n) - 



(E17) 
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We have explicitly separated the medium specific contributions from the vacuum 
projections. Medium specific contributions are either condensates that contain the 
medium four-velocity v or combinations of condensates that appear in vacuum and 
medium. The latter ones are always written with angled brackets. Applying vacuum 
projections to the medium specific terms makes them zero in the vacuum limit. 

In order to give an example for this procedure we briefly proof Eq. (1E16|) . Due 
to Lorentz covariance we write 

{^\ql5laQ^.uq\^) = At^^.^v'' . (E18) 
Our aim is to determine the Lorentz scalar A. We can write for any four-vector 

va^'G^u = ^t;a7"7V6^M- (E19) 
Expanding the product of Dirac matrices in terms of the Clifford algebra, one obtains 

7"7^" = 9""^ + g'^^r + g'^^l" + ^e™'^"757<. • (E20) 

Due to the equations of motion and by the definition of the gluon field strength 
tensor, one can show that 

{n\dv^'-i''g^^d\n) = o. (E2i) 

Altogether, this gives the relation 

{n\dva^''g,,d\n) = -{n\dj,j^g,,d\n)e''''''''v^. (E22) 



Contracting (1E18|) with e"A*^<^ and using e^^'^'^e = Gg'''^ , one can show the desired 



relation flElGl) . 

Let us define the operator T in such a way that it creates a traceless expression 
with respect to Lorentz indices. For two Lorentz indices T reads 

fiO,,) = 0,,-^02, (E23) 



and we immediately observe that the medium specific terms in ( lEllI) and in (IE14P 



originate from the contraction of v^v'^/v'^ with a certain traceless expression: 
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{n\qf (7^:^,) q\n) 



{n\qf q\n) 



4 

777,2 -1 

-^{Q\qq\Q) - -{Q\qg^agq\n) 



(E24) 



+{n\q^-Lq\n) 



(E25) 



This is clear since the additional medium contributions can be obtained by per- 
forming the complete projection of the structure to be projected minus the vacuum 
projection, which is proportional to products of the metric tensor, O^j, — ^g°'^Oap- 
Therefore, the additional medium contributions must be traceless. The generaliza- 
tion to arbitrary Lorentz indices is obvious. 
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